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ABSTRACT 

Star formation generally proceeds inside-out, with overdense regions inside protostellar cores col- 
lapsing rapidly and progressively less dense regions following later. Consequently, a small protostar 
will form early in the evolution of a core, and collapsing material will fall to the protostellar surface 
and radiate away its gravitational potential energy. The resulting accretion luminosity will heat the 
core and may substantially affect the process of fragmentation. This is of particular interest for mas- 
sive cores that, at their initial temperatures, have masses much greater than a thermal Jeans mass 
and thus might be expected to fragment into many stars during collapse. Here I show that accretion 
luminosity can heat the inner parts of a core to > 100 K very early in the star formation process, 
and that this in turn strongly suppresses fragmentation. This has implications for a number of out- 
standing problems in star formation, including the mechanism of massive star formation, the origin 
of the stellar initial mass function and its relationship to the core mass function, the demographics of 
massive binaries, and the equation of state in star-forming gas. 

Subject headings: binaries: general — equation of state — ISM: clouds — methods: numerical — 
radiative transfer — stars: formation 



1. INTRODUCTION 

In the last few years, millimeter and sub-mm obser- 
vations have reached resolutions sufficient to identify 
dense molecular condensations with masses ~ 100 AIq 
and sizes ^ 0.1 pc in regions o f massive sta r forma- 
tion (a recent review is given in IGaravl l200l . These 
objects are cold, supersonically turbulent, and centrally 
condei ised. In some cases they show signs of coherent 
infall ijBeuther et alJl2005|) . The condensations have a 
mass distributi on that resembles the stellar initial mass 
function (IMF, iBeiither fc Schilkel l20(M': 'Rcid & Wilson" 
I200 5'). which naturally leads to the question of whether 
they could collapse to form massive stars, making them 
the high mass analogs of the cores from which low mass 
stars form. 

One objection to this scenario is the possibility of 
fragmentation. Since the thermal Jeans mass in a cold 
(T - 10 K), dense {uh 10^ cm'^) gas is - 1 M©, 
why should a core of mass ~ 100 Mq collapse to form 
a single massive object object (or a sm all-multiple sys- 
tem), rather than many smaller objects? iBate & BonnelM 
argue based on simulations that denser regions 
produce smaller fragments, since the Jeans mass de- 
creases with density at fixed tem perat ure (although see 
iMartel et alJ200fij^ . iDobbs et all pOO^D simulate the col- 
lapse of a massive turbulent core and find that it frag- 
ments to as many as 20 and as few as 2 objects, depend- 
ing on the assumed equation of state. If fragmentation 
of massive cores to many objects were common, there 
would be no direct mapping between the core and stellar 
mass distributions; their agreement would simply be a co- 
incidence. If, on the other hand, massive cores typically 
collapse to one or a few objects, then one can explain the 
IMF in terms of the core mass function. (Fragmentation 
to 2 — 3 objects does not pose a problem in mapping 
from core to star masses, because the massive star IMF 
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is uncorrected for multiplicity.) 

Almost all published work on fragmentation of massive 
cores treats the gas as either isothermal or barotropic. 
The barotropic approximation is based on calculations 
showing that radiative cooling keeps gas isothermal until 
it reaches a density of ~ 10"^"^ g cm~^ at which point 
it tra i isitions to adiabatic llLarsonI 1969HMasunaga et alJ 



119981 iMasunaga fc Inutsukal l200m . However, this ap- 
proach assumes that gas is heated solely by compression. 
While this is probably reasonable far from point sources 
of radiation, it is likely to fail once a protostar forms, 
since heating of the gas by irradiation from the lumi- 
nous, accreting central object vastly exceeds heating due 
to gas compression. Indeed, in detailed one-dimensional 
radiation-h ydrodynamic calculations of t he evolution of a 
1 Mq core. IMasunaga fc Inutsiikal l)2000t) find that accre- 
tion lu minosity can heat gas to > 100 K out to hundreds 
of AU.'Matzn er fc Leviiil l)2005j) find that heating due to 
accretion luminosity is sufficient to prevent fragmenta- 
tion of low mass protostellar disks into brown dwarfs. 

Massive star-forming regions have high densities that 
produce accretion rates orders of magnitude larger than 
in low mass regions, and they are optically thick even 
in the mid-infrared. Thus, heating by accretion lumi- 
nosity within them is likely to have even more profound 
effects. In this paper I test the assumption that mas- 
sive cores can be described with a barotropic equation of 
state, and explore how a more realistic treatment of gas 
thermodynamics changes the picture of fragmentation in 
these objects. In § [3 I develop a simple analytic model 
for the temperature distribution in massive cores, and in 
§ 121 1 use the results to estimate how fragmentation is 
likely to be affected. Finally, in §01 1 discuss some of the 
broader implications of my findings. 

2. THE TEMPERATURE STRUCTURE IN MASSIVE CORES 

Since the focus of this paper is fragmentation, I concen- 
trate on cores at early times when the gas mass greatly 
exceeds the stellar mass, and accretion is the dominant 
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source of luminosity. To make tlie calcula tions definite, I 
use th e turbulent core model proposed bv lMcKee fc TanI 
l)2003l hereafter MT03) , in which the density distribution 
within a core is a power law, p cx r^'^p , and both MT03 
and I adopt fcp = 1.5 as a fiducial value. In this model 
a choice of core mass rric and mean column density Sc 
specifies the outer radius of the core = [m^/ {-kYjc)]^^"^ 
and hence the density profile. Requiring hydrostatic bal- 
ance gives the effective sound speed versus radius, a = 
ao{r/roY-''^'\ with = [2{kp - l)]-^/\GmJrof/\ 
Within the core is an embedded protostar of mass m,. 
This model neglects turbulent structure within a core, 
the effects of which I discuss in more detail in § O 

One can estimate the accretion rate in a massive core 
from equation (41) of MT03: 

(1) 

(This calculation assumes that half the mass reaching the 
star is blown away in an outflow.) For a protostar in the 
mass range 0.01 — 1 Mq, this gives a typical accretion rate 
of roughly 10~^ — 10"'* Mq yr~* for a core with TOc = 50 
Mq and Sc = 1.0 g cm~^. Though I compute m* from 
MT03, one can obtain similar values by simple order- 
of-magnitude arguments. For example, iBeuther et al.l 
((20051) observe a velocity dispersion a = 1.0 km s~* in 
IRDC 18223-3 main, giving an e xpected accretion rate 
~ a^/G « 10-4 Mq yr~\ \Reid fc Wilsm] ll200l 
find cores with typical values rric ^ 100 Mq, To ~ 0.25 
pc. At the time when the forming protostar has mass m, , 
the accretion rate should be rh* ~ TO*/tff(m*) ~ 10~^ 
Mq yr~*, where tff(m*) is the free-fall time at the edge 
of the central region in the core that contains m, of gas, 
and the evaluation assu mes kn = 1-5 . Simulations ob- 
tain comparable values. tPobbs et al. (2005) find accre- 
tion rates onto single objects of 10~^ — lO"** Mq yr^^ 
for protostars of mass m, < Mq in a core with rric = 30 
Mq, = 0.6 g cm-2. 

To compute the accretion luminosity, one must know 
the protostellar radius. Initially collapsing gas forms a 
pressure-supported "first core" a few AU in size, but this 
becomes unstable and collapses to stellar sizes once its 
mass reaches roughly 0.05 (jMasunaea et alJ Il998t 
IMasunaea fc Inuts uka| I200(]|) . Thereafter one can com- 
pute the protostellar radius using the simple evolution 
model of MT03, which agrees well with the more de - 
tailed numerical calculations of IPalla fc Stable^ l(1992j) . 
The accretion luminosity is Lace = /accWH.77i*/rH., where 
/acc ~ 0.5 is the fraction of the accretion energy that goes 
into radiation rather than driving a wind. For typical 
massive star-forming regions, this gives iacc 10 — 100 

Lq. 

Given a central luminosity and a density structure, one 
could numerically determine the temperature structure 
(e.g. I^vczic & Elitzur 1997; Whitncv et al. 2003). How- 
ever, for simplicity I in stead use the analyt ic approxima- 
tion given bv tChakrabarti fc McKeel l)2005() . which agrees 
well with numerical calculations. I will not review the full 
Chakrabarti & McKee formalism, but the central idea is 
that the majority of the energy radiated away at a given 
frequency comes from a thin radial shell whose position 
is determined by a competition between opacity, which 
suppresses emission from small radii, and decreasing tern- 
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Fig. 1. — Temperature versus radius for models with the val- 
ues of rric and Sc indicated in the panels. Temperatures shown 
are for protostellar masses of m« = 0.8 Mq, m, = 0.2 A/q, and 
rrit = 0.05 Mq (solid lines, top to bottom) and for the Dobbs et al] 
120051) barotropic equation of state (dashed line). I also show 
the effective temperature set by the turbulent velocity dispersion, 
Tturb = 2.33 m^f CT^/fc^ (dotted line). 

perature, which reduces emission from large radii. Once 
one knows the characteristic emission radius versus fre- 
quency, one can estimate the temperature profile, which 
is roughly a power law T = Tch{r / Rch)^'^'^ ■ For a given 
model, I compute Rch, Teh, and kx from Chakrabarti & 
McKee 's equatio ns (6), (7), an d (41 ), using the opacity 
law of iWeingartner fc Drain3 l|2001f) . Since I am con- 
cerned with early times when the protostar has accreted 
a negligible fraction of the cloud, I assume that the core 
density profile is unchanged from its initial MT03 form. 

I plot the temperature profile for a selection of pa- 
rameters in Figure ^ The gas is clearly quite warm, 
but the true tempera t ure is likely even hotter. The 
iChakrabarti fc McKed l|2005f) approximation underesti- 
mates the temperature at very small radii, I have ne- 
glected compressional heating, and I have assumed that 
all of the energy that goes into driving a wind escapes 
from the core rather than being radiated away by shocks 
within it. Figure ^ also shows the te mperature set by th e 
barotropic equation of state used bv lDobbs et al.l l|2005f) . 
T = 20 K for p < 10-^'^ g cm'^, T = 431 K for p > lO'^^ 
g cm~'^, T oc p'^^^ in between , which is siniilar to ones 
used in other simulat i ons (e.g.lBonnell et alJl2003l 120041: 
IBonnell fc Batai2005HBate fc Bonnel]ll2005D . As the fig- 
ure shows, the barotropic approximation severely under- 
estimates the temperature. One can understand this 
intuitively: the barotropic approximation says that, at 
some critical density, gas becomes adiabatic and ceases 
radiating away the energy it gains through compression. 
However, since potential energy varies as r~^, the vast 
majority of the energy released by a collapse comes out 
in the final plunge onto the stellar surface, and is then 
radiatively transferred to the rest of the gas. Even when 
the mass accreted onto a star is a small fraction of the 
core mass, the energy released by its fall to very small 
radii can dominate the total gravitational potential en- 
ergy released. By including only the energy released by 
the fall from to ^ 10 AU, roughly the numerical res- 
olution of most simulations, and not the energy released 
from '^10 AU to ~ Rq, the barotropic approximation 
misses the vast majority of the energy of collapse. 

To validate my use of the equilibrium temperature pro- 
file, I also estimate the time required to reach radiative 
equilibrium. If -Ediff' is the difference in gas thermal en- 
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ergy between a core whose temperature is as shown in 
Figure Q a-nd a core with a constant temperature 20 
K, the time to reach radiative equihbrium is roughly 
irad = E^iff/Liicc- This estimate neglects radiative cool- 
ing, but since the cooling rate varies as T"* in optically 
thick gas, cooling is negligible compared to heating until 
the gas is very close to equilibrium. For the parameter 
range = 0.05 - 1.0 Mq, = 50 - 150 Mq, and 
Tjc = 0.3 — 1.0 g cm^^, iiad < 10 yr. This is shorter than 
the free-fall time everywhere except in the central few 
tenths of an AU, so the assumption of radiative equilib- 
rium is justified. 

3. FRAGMENTATION IN MASSIVE CORES 

The high temperatures produced by accretion luminos- 
ity will change how fragmentation proceeds. A cold cen- 
trally condensed core will undergo global collapse, but 
local inhomogeneitics produced by turbulence will also 
be gravitationally amplified, producing a collaps ing clus- 
ter of small stars, as seen in lDobbs et al.l ()2005j) . rather 
than monolithic collapse to a single star. Heating pre- 
vents the growth of small perturbations, thereby both 
reducing the number of perturbations that grow and in- 
creasing the mean mass of those that can. To study this 
effect, I compute the Jeans length Aj — [■Kcl/{Gp)]^^'^ 
as a function of radius r in my model cores. (I take 
the mean particle mass to be 2.33m//.) At radii where 
Aj > 2r, no perturbations that fit within the sphere 
of radius r are Jeans unstable and no fragmentation 
can occur. The radius where Xj — 2r defines a length 
scale rtiag and mass scale mfrag below which there can 
be no fragmentation. I compute the Jeans length using 
the thermal sound speed rather than the effective sound 
speed set by turbulence because both simulations and 
theory show that the typical fragment mass in a tur- 
bulent medium is the thermal Jeans mass rather than 
the "turbulent" Jean s mass (Padoan & Nordlund 2002; 
iMac Low &: Klessenl2004l and references therein) . Phys- 
ically, this occurs because fragment formation in a tur- 
bulent medium tends to occur at stagnation points that 
are relatively non-turbulent. 

Figure 121 shows rfrag and mfi-ag for a range of parameter 
choices. For a barotropic equation of state, there is no 
heating, so the mass scale is constant. For rrir = 30 M(t,. 
S = .6 g cm^-^, the parameters used bv iDobbs et al.l 
l)2005|) . TOfrag = 0.3 Mq, which roughly agrees with the 
typical fragment mass in the simulations. However, when 
one includes radiation, the resulting high temperatures 
suppress fragmentation to objects of mass ^ 1 Mq over 
length scales ^ 1000 AU from the moment that the first 
core collapses and a protostar appears. The barotropic 
approximation underestimates fragment masses by fac- 
tors as large as 10, and it switches the problem from 
the regime where Wfrag > (which favors monolithic 
collapse) to Wfrag < 1^* (which favors fragmentation). 

As illustrated in Figure^ the velocity dispersion in the 
radiatively heated gas within a few hundred AU of the 
star can be larger than the velocity dispersion required to 
maintain hydrostatic balance. As a result, the gas may 
temporarily expand, until infall from the colder outer 
parts of the core raises the density. This could further 
suppress fragmentation. If the expansion is adiabatic 
with 7 = 5/3, the Jeans length will increase as A,/ cx 




Fig. 2. — Mass mf^ag and radius rfj-ag below which fragmentation 
is suppressed for (mc,Sc) = (5OM0,O.3 g cm~^) {thin lines) and 
(150 M0, 1.0 g cm~^) (thick lines). I plot mfi-ag {upper left), rf^ag 
{lower left), the ratio of mf^^g computed using radiative models 
to the value one infers using a barotropic equation of state (upper 
right), and the ratio of mf^ag to the protostellar mass m* (lower 
right). For all plots but the ratio of mf^ag computed radiatively 
and barotropically, I show both radiative transfer (solid lines) and 
barotropic (dashed lines) models. The dotted vertical line at r?i» = 
0.05 Mq indicates the mass at which a protostar first forms. 

With radiative effects, the expansion effectively 
decreases Ec at constant nic, the net effect of which is 
also to increase Aj. 

Finally, note that turbulence in massive cores likely 
produces multiple density peaks within the overall 
centrally-condensed structure. The densest of these 
peaks (which is almost certainly close to the center) has 
the shortest free-fall time and will collapse first. At the 
accretion rate given by equation for rric = 50—150 
Mq and Ec = 0.3 - 1.0 g cm"^^ o.05 Mq of gas wiU accu- 
mulate, collapse to stellar density, and begin irradiating 
its surroundings in 30 — 300 yr, which is vastly shorter 
than the mean-density free-fall time of 3 x 10* — 1 x 10^ 
yr. Unless there is an unreasonably synchronized collapse 
of several density peaks, all the gas in the cloud except 
the first 0.05 Mq to collapse will be subject to radiative 
feedback. In a turbulent core self-shielding may leave 
sufficiently overdense regions cooler than the spherical 
calculations indicate, so fragmentation can only be fol- 
lowed in detail by radiation-hydrodynamic simulations. 
Nonetheless, the spherical calculations definitively show 
that radiative heating will raise the Jeans length and 
mass over a large volume within the core well before the 
vast majority of the gas collapses, so that the number of 
fragments should be significantly fewer, and their mean 
mass much larger, than one would estimate by assuming 
an isothermal or barotropic equation of state. 

4. IMPLICATIONS AND CONCLUSIONS 

4.1. Massive Stars, Cores, and the IMF 

The realization that radiative heating can suppress 
fragmentation argues that massive condensations really 
are cores, in that they are likely to collapse to pro- 
duce single stars or small multiple systems, not clus- 
ters. This supports the idea that the massive stars 
form by accretion (MT03) and that the mass func- 
tion of stars is_ determined at the stage of fragmenta- 
tion into cores fPadoan & Nordlund 2002; Beuthcr et al] 
120051 IReid fc W ilson 2005). It removes the need to ex- 
plain the common shape of the core and stellar mass 
functions as pure coincidence. 
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In contrast, the suppression of fragmentation by ra- 
diative heating presents a ser ious problem fo r mod- 
els of the type proposed by 'Bat e fc Bonnelll l)2005f) 
and iBonnell fc B ate (2005) in which ejection and com- 
petitive accretion determine the IMF and direct coUi- 
sions are required to create massive stars (though see 
iKrumholz et al.ll2005l for a critique of these models on 
other grounds). In such models, all stars are born as 
brown dwarfs or very low mass stars in clusters ~ 1000 
AU in size. Most objects then accrete until being ejected 
by N-body interactions, while a few experience run- 
away collisions, leading to the formation of massive stars. 
However, in dense regions where the accretion luminosity 
is high, the formation of the first protostar within a core 
will inhibit the formation of subsequent ones around it. 
Rather than producing a large cluster of low mass proto- 
stars or brown dwarfs as the competitive accretion and 
collision models demand, heating will produce a small 
number of larger protostars. 

It is unclear whether these results present a problem 
to IMF models that depend on the structure of the equa- 
tion of state for collapsing gas l)Larsonll2005|) . If the typ- 
ical separation between the fragments in a protocluster 
is larger than ~ 1000 AU, then, as Larson suggests, ra- 
diative heating will not substantially modify the typical 
fragment mass until enough stars have formed to heat 
the entire protocluster. On the other hand, if the typical 
inter-fragment spacing is small, then radiative feedback 
effects cannot be ignored. 

4.2. Binary Formation 

Suppression of fragmentation by feedback also provides 
a natural explanation for the observation that, when 
massive stars are in binaries, the binaries tend to consist 
of two massiv e stars rather than a high mass star and a 
low mass one l|Pinsonneault fc Stanek 2006^1 . If binaries 
result from the fragmentation of massive cores, then one 
would expect the resulting stars to be massive because, 
once the first protostar forms, it will raise the Jeans mass 
and thereby favor the formation of massive stars there- 
after. Although I have only considered heating due to ac- 
cretion luminosity, it is likely to be even more significant 
once protostars become massive enough to start fusing 
deuterium and then hydrogen. This will greatly increase 



their luminosity, and raise the Jeans mass even more. 
Testing whether this effect can quantitatively explain 
the observed properties of binaries will require radiation- 
hydrodynamic simulations. Radiation f eedback ma y also 
be important for low mass binaries. iLadal |2006) re- 
cently pointed out that most low mass stars are single- 
tons rather than binaries. Barotropic simulations of low 
mass star formatio n almost inv ariably produce multiple 
systems ( Goodwin et al.lf2004|) . and radiation feedback 
may provide a way of reducing the multiplicity. 

4.3. The Barotropic Approximation in Simulations 

In some cases the barotropic approximation produces 
results at odds with radiative transfer calculation s even 
when there is no point source of radiation (Bos s et alJ 
l200fl() . Feedback makes the problem far worse. Simula- 
tions on scales 1000 AU, such as those that model core 
formation within a larger molecular c loud or clump (e.g. 
iLi et a l. 2004: Tillcv & Pudritz 2004), arc probably safe. 
On the other hand, as Matzncr & Levin (2005) point out 
for disk fragmentation, the barotropic approximation is 
likely to produce incor rect results in s imulations t hat go 
to sm aller scales (e.g. iBate et all 1200 3: Goodwi n et alJ 
'2004^. For this reason, simulations of fragmentation 
should be done with radiative transfer. However, in sim- 
ulations that do use a barotropic equation of state, one 
can at least attempt to check the validity of the approx- 
imation ex post facto by selecting a few time slices and 
calculating the temperature distribution, including ac- 
cretion lumino sity, using a radiative transfer code (e.g. 
IWhitnev et aLl li2003). If the temperature is significantly 
higher than indicated by the assumed equation of state 
in the simulation, then the results for fragmentation in 
the region whose temperature has been underestimated 
must be regarded with suspicion. 
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